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Abstract 

We present a new kind of particle path in constant vorticity water of finite depth, 
within the framework of small-amplitude waves. 

1 Introduction 

A peakon is a soliton with discontinuous first derivative |17j . The concept was 
introduced in 1993 by Camassa and Holm in the paper [4], where they derived 
the CH shallow water equation 

Ut + 2KU X + 3UU X — Utxx = 2u x Uxx + uu XX x, (CH) (1) 

(x,t) e R x (0, oo), K being a real constant. Alternative derivations of CH 
equation are provided in the papers [22] , [9] , [19] . The peakons arise as solution 
of this equation for n=0. The CH peakons are given by 

u(x,t) = cexp(— | a; — ct\), c € R. (2) 

Since peakon solutions are only piecewise differentiable, they must be inter- 
preted in a suitable weak sense. The derivative 

u x = — c sgn [x — ct) exp(— \x — ct\) (3) 

has a jump discontinuity at the peak. The second derivative u xx must be taken 
in the sense of distributions and will contain a Dirac delta function 

u xx — cexp(— \x — ct\) — 2cS(x — ct) (4) 

The function m is defined by 

m(x, t) := u — u xx = 2c5(x — ct). (5) 
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Physically m has the interpretation of momentum [?] , [TB] . 

The peakon ([2]) has amplitude c and travels at speed c. At x = ct the momentum 

§5$ blows up at +oo. 

A small perturbation of a CH peakon yields another one which remains close 
to some translate of the initial one at all later times. In this sense the CH 
peakons are orbitally stable [TU]- Of particular interest is the description of 
peakon dynamics in terms of a system of completely integrable Hamiltonian 
equations for the locations of the peaks of the solution. Thus, each peakon 
solution can be associated with a mechanical system of moving particles. Being 
solitons, they retain their shape and speed after interacting with other peakons 
PQ. The peakon interaction plays an important role in the general dynamics 
of the solutions to the equation (see the discussion in [16]) and provided the 
framework for the construction of global weak solutions both in the conservative 
case [5] as well as in the dissipative case [3] . One of the main interests in CH 
equation was that, in contrast to other standard shallow water equations, as for 
example the KdV equation, it models breaking waves: smooth solutions that 
develop singularities in finite time, the solution being bounded but its slope 
becoming unbounded. This fact was already noted in [4] and subsequently 
proved in [BJ. 

Another completely integrable CH-type equation which has peakon solutions 
P2] is the Degasperis-Procesi equation [TJ] 

u t + Auu x - u txx = 3u x ii xx + uu xxx , (DP) (6) 

(x, i)eRx (0, oo). The DP equation possesses not only peaked solitons ([2]) but 
also discontinuous solitons, so-called shock-peakons [25] of the form 

u(x,t) = cexp(-|a;|) - — !— sgn (x) exp(-|x|), k>0 (7) 
t + k 

At the peak they have a finite jump in the function u itself. The shock-peakon 
solutions must be interpreted in a proper weak formulation. The derivative u x 
will contain 8 and the function m := u — u xx will be a linear combination of 8 
and 8' distributions. It is not known to the author if the function m can be in 
this case interpreted as momentum. We point out that the CH equation with 
k = 0, K 7^ 0, is a geodesic equation on the diffcomorphism group of the circle 
[8], respectively on the Bott-Virasoro group [26], [7], while the DP equation is 
a non- metric equation |15j . 

The shock-peakon moves at constant speed c (in particular, does not move 
if c = [25]) which is equal to the average amplitude at the jump. The shock 
"dissipates away" like 1/t as t — > +oo. 

The peakons of the DP equation are also true solitons that interact via elastic 
collisions under the DP dynamics [23], and are also orbitally stable [23 . 

In what follows we will see that in the study of particle motion beneath small- 
amplitude water waves in constant vorticity flows a peakon trajectory comes up. 
This solution contains arctanh(-) function, having a vertical asymptote in the 
positive direction. 
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2 Particle path beneath small-amplitude water waves in 
constant vorticity flows 

We consider two-dimensional gravity waves on constant vorticity water of finite 
depth. They are described, in non-dimensional scaled variables, by the following 
boundary value problem (see, for example [21j): 

u t + e(uu x + vu z ) = -p x 
S 2 [v t + e(uv x + vv z )} = -p z 
u x + v z = 

u z = S 2 v x + ^uj n (8) 
v = rf t + eurj x on z = 1 + er](x, t) 

p = i] on z = 1 + erj(x, t) 

v = on z = 

where (x, z) are the space coordinates, (u(x, z,t),v(x, z, t)) is the velocity field 
of the water, p(x, z,t) denotes the pressure, g is the constant gravitational ac- 
celeration in the negative z direction, ujq being the constant vorticity. We have 
introduced the amplitude parameter e = and the shallowness parameter 

5 = ^j-, with a the amplitude of the wave and A the wavelength, ho > is 
the undisturbed depth of the fluid and z — 1 + er)(x, t) represent the free upper 
surface of the fluid in non-dimensional scaled variables. The existence of solu- 
tions of large and small amplitude was recently proved in [11] where it is also 
shown that linearization provides an accurate approximation for waves of small 
amplitude. 

By letting e — > 0, 5 being fixed, we obtain a linear approximation of our 
problem, that is, 

u t + Px = 
S 2 v t + Pz =0 
u x + v z = 

u z = S 2 v x + ^f^ (9) 
v — rjt on z — 1 

p = rj on z — 1 

v = onz = 

The system (0) has the solution 
rj(x. t) — cos(2-7r(a; — ct)) 

< x > z > *) = ^ts) cosh(27rrJz) 003(2^(2: - ct)) + ^SK Z + Co 

z ' *) = sinh^i) sinh ( 27rfc ) sin(27r(x - ct)) 
P( x > z ' *) = Bi nh(2^) cosh(27T(5z) cos(27r(a; - ct)) 

with the non-dimensional speed of the linear wave given by 

9 tanh(27r<5) 



(10) 
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Let (x(t),z(t)) be the path of a particle in the fluid domain, with location 
(x(0), z(0)) := (xo,z ) at time t — 0. The motion of the particles below the 
small-amplitude gravity water waves given by (fit)]), is described by the following 
system of differential equations 



f = u(x, z, t) = cosh(27n5z) C0S (2ir(x - ct)) + 



CO 



f t = V ( x , Z} t) = sinh(2^z) sin(2^(x - ct)) 

(12) 

Notice that the constant cq is the average of the horizontal fluid velocity on the 
bottom over any horizontal segment of length 1, that is, 



Co 



x + l 



u(s, 0, t)ds 



(13) 



This is accordance with Stokes' definition of the wave speed for irrotational flows 
(see the discussion in |12jh 

To study the exact solution of the system (fT2")) it is more convenient to rewrite 
it in the following moving frame 



X = 2n(x - ct), Z = 2ir5z 
This transformation yields 



dX 


4tt 2 Sc 


dt 


sinh(27T(5) 


dZ _ 


47r 2 <5c 


dt 


smh(27r<5) 



We denote by 



A := 



47T 2 (5c 



and CI 



sinh(27T(5) gS 

With the notations (fTB)) , the system (|T5|) becomes: 



dX 
dt 



Acosh(Z) cos(X) + ft Z + 27t(c - c) 



^ = A sinh(Z) sin(X) 
We write the second equation of this system in the form 

dZ 



sinh(Z) 



A sin X{t)dt 



Integrating, we get 



log 



tanh 



J AsmX(t)dt 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 
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If 



then 



J AsinX(t)dt < 



exp ^ AsinX(t) dt 



Z(t) = 2 arctanh 
Taking into account the formula: 

1 + tanh 2 (x) 



cosh(2;r) = 



(20) 
(21) 

(22) 



1 — tanh 2 (x) ' 

and the expression (|21|) of Z(t), the first equation of the system (|17|) becomes 



dX 1 + w 2 

—— = A- ~ cos(X) + 2f2o arctanh (w) + 2tt(co — c) 

dt 1 — w z 



where we have denoted by 



■;(*) := exp (J A sin X(t)dtj 



With (|20|) in view, we have 
From ([2l]) we get 



< w < 1 



AsinX(t) = 



1 dw 



w(t) dt 

Differentiating with respect to t this relation, we obtain 



,dX 1 
A cos X — = — 
dt w 2 - 



dt 2 



■w 



dw 
~dt 



From (1261) we have furthermore 



1 / dw 



A z cos z (X) = A 2 , . 

w z \ dt 

Thus, taking into account (|2"7|) . (|2"5|). the equation (f2"3"|) becomes 



(23) 

(24) 

(25) 
(26) 

(27) 
(28) 



2w / dw\ 



dt 2 l~w 2 \dt J 



1 + w 2 

I 

1 — w 2 



2fi arctanh (iu) + 2tt(c - c) =0 (29) 



We make the following substitution 

( 2 (w) :=A 2 w 2 



dw\' 
~dt 



(30) 
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Differentiating with respect to t this relation, we get 

tUL = A 2 w - — 
dw dt 2 

We replace f|30[) . ([31]) into the equation (|29l) and we obtain the equation 



2w 



dw 1 



-e 2 



2f^o arctanh (w) + 27t(cq — c) 



A solution of the equation (|3"2l is 

e = o 

which, in view of ([30]) and (|26|) implies 

sinX(i) = ±1 



(31) 

(32) 

(33) 
(34) 



Therefore, from ([2~Tj) with the condition (I2TJ1) . and further from (fT4^1 . a solution 
of the system (fT2"]) is 



a;(i) = c£ + fci 

^(t) = arctanh [exp (— A\t\)] 

k\ being a constant. We observe that 

lim x(t) — fci, lim z(t) = lim z(t) = +oo 
t-tO t->0 t-yO 

t>0 t<0 



and 



lim = ±oo, lim z(t) = 



(35) 



(36) 



(37) 



Therefore, a; = k\ will be a vertical asymptote and z — will be a horizontal 
asymptote for the curve (|35|) . The graph of the parametric curve (|35[) is drawn 
in figure 1. 



Figure 1. 




Notice that within the setting of irrotational flows with no underlying current 
(see [5]) there are no such paths but in the context of irrotational flows with an 
underlying (uniform) current, the possibility of such paths was already noticed: 
see [12] for the exact solutions, where this shape can be thought of as a limiting 
case of the situation depicted in Figure 4.4 (ii), as well as [20] for the linearized 
problem, where somewhat similar particle paths are encountered. 
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